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We report the observation of strongly temperature-dependent, asymmetric spectral lines in elec- 
tronic Raman scattering of graphite in a high magnetic field up to 45 T applied along the c-axis. The 
magnetic field quantizes the in-plane motion, while the out-of-plane motion remains free, effectively 
reducing the system dimension from three to one. Optically created electron-hole pairs interact 
with, or shake up, the one-dimensional Fermi sea in the lowest Landau subbands. Based on the 
Tomonaga-Luttinger liquid theory, we show that interaction effects modify the van Hove singularity 
to the form {to — A) "~ ' at zero temperature. We predict a thermal broadening factor that in- 
creases linearly with the temperature. Our model reproduces the observed temperature-dependent 
line-shape, determining a to be ~0.05 at 40 T. 

PACS numbers: 78.30.-j, 71.70.Di, 73.61. Cw, 76.40.+b, 78.20.Bh 



Electron-electron interactions are progressively more 
important as the system dimension is lowered. One- 
dimensional (Id) systems, in particular, provide model 
environments in which to explore interaction effects [l|. 
Interacting Id electrons are expected to form an exotic 
electronic state of matter, the Tomonaga-Luttinger liq- 
uid (TLL) [245|- A strong magnetic field, B, can sup- 
press the electrons' kinetic energy, thus enhancing the 
effect of interactions, as exemplified by the fractional 
quantum Hall effect [6|-|8|. In a 3d material, an applied 
magnetic field creates an effective Id system along the 
field, ideal for a systematic study of interaction effects in 
a highly controllable fashion [9|. Particularly promising 
are 3d metals with small electron and/or hole pockets 
near the Fermi energy (Ep), such as bismuth |10l4l4| and 
graphite 1^, llSMlSfl , where the magnetic quantum limit 
can be readily reached with B '--^lO T. 

Here we use Raman spectroscopy to study electronic 
states and correlations in natural graphite in a strong 
B up to 45 T, applied along the c-axis. This quantizes 
the electronic motion in the afc-plane, while the motion 
along the c-axis remains free, thus reducing the effec- 
tive dimension from three to one. Instead of the main 
Raman features related to phonons [19l |20| , here we fo- 
cus on a series of electronic Raman features assigned to 
electronic inter-Landau- level (LL) transitions [2l|, whose 
B-dependence can be explained thr oug h the Slonczewski- 
Weiss-McClure (SWM) model [2i - [24 |. Each feature ex- 
hibits strongly temperature (r)-dependent shape. Our 
calculations show that scattering by thermally excited 
acoustic phonons 2514281 is too weak to explain the ob- 
servations. Electron-electron interactions, on the other 
hand, are shown to be the cause for the observed T depen- 



dence, through the 'shake-up' process known in the prob- 
lem of X-ray (or Fermi-edge) singularities [5|. Namely, 
optically created electron-hole pairs interact with, or 
shake up, the Id Fermi sea in the lowest Landau sub- 
bands, resulting in line-shape deviations from single- 
particle densities of states (i.e.. Id van Hove singular- 
ities). Based on the TLL theory [lH5|, we show that 
electron-electron interactions modify the van Hove sin- 
gularity to the form {uj — A)^"~^'^ at K, where uj is 
the photon frequency, A the band-edge frequency, and 
a a dimensionless measure of the influence of electron- 
electron interactions. At finite T, we predict a thermal 
broadening factor, T{T) ex T. Our model reproduces the 
observed T-dependent line-shape, determining a to be 
0.016, 0.026, and 0.048, at 20, 30 and 40 T, respectively. 

Raman measurements are performed on natural 
graphite (NGS Naturgraphit GmbH) in a back-scattering 
Faraday geometry in B up to 45 T, as described in 
Ref. |2lJ. The excitation beam from a 532-nm laser is fo- 
cused to a spot size of <20 //m with a power of ~13 mW. 
Most of the data are collected with a spectral resolution 
of -3.4 cm-i. For the high-B, low-T (<10 K) mea- 
surements of the sharpest peaks, a spectral resolution 
of —1.9 cm~^ is used. The temperature drift over an 
integration time of up to 7 minutes, monitored by a tem- 
perature sensor installed below the sample, is <1 K at T 
= 7 K and <2 K at T > 180 K. 

Figure [TJa) shows Raman spectra taken at 10, 20, 
and 30 T at 7 K. The main feature is the G peak at 
—1580 cm~^, due to E2g phonons [20|. In the presence of 
B, electronic Raman features appear, coming from inter- 
LL transitions, labeled (1,1), (2,2), ... etc., which we 
focus on in this work. Figure [TJb) shows a series of spec- 
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FIG. 1: (Color online) (a) Raman spectra at 10, 20, and 
30 T. The feature at ~1580 cm"^ is the G-peak due to the 
E2g phonons. Inset: schematic energy level diagram showing 
the transitions responsible for the electronic peaks, (b) Data 
taken at various B at T = 7 K, showing peaks due to (1,1) 
through (4,4) interband transitions, (c) Peak positions of the 
observed (1,1), (2,2), and (3,3) transitions as a function of B, 
together with calculations based on the SWM model. 



tra taken at various 5 at 7 K, exhibiting electronic peaks 
that move with B. These peaks can be attributed to the 
"symmetric" inter-LL excitations in the vicinity of the K 
point [2l|,|29|. The strongest, lowest-frequency transition, 
is (1,1), which is from the n = — 1 level in the valence 
band to the n = I level in the conduction band. Similarly, 
we observe the (2,2), (3,3), and (4,4) transitions [see also 
the zero-field in-plane dispersions and energy levels near 
the K-point in the inset to Fig. [Ua)]. The symmetric 
inter-LL excitations are non-resonant Raman processes 
and have been theoretically investigated for single-layer 
graphene (SLG) [31| and bilayer graphene (BLG) |33 |. 
The peak positions of the three lowest-energy transitions 
are plotted against B in Fig. [TJc); they agree well with 
our calculations (solid and dashed lines) based on the 
SWM model ^. 

These inter-LL transitions show strong T dependence, 
as indicated in Fig. O where Raman spectra at various 
T are plotted for (a) 20, (b) 30, and (c) 40 T. At the 
lowest T, the peaks exhibit sharp and asymmetric line- 
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FIG. 2: (Color online) T-dependent electronic Raman scatter- 
ing of graphite at (a)20 (b) 30, and (c)40T. (d) T dependence 
of the broadening factor, P, of the (1,1) line at 20 T (solid 
circles) and 30 T (open circles). The lines are fits to the data 



shapes, reminiscent of a Id van Hove singularity, as ex- 
pected from the effective dimension reduction from three 
to one in a B. As T increases, there is significant peak 
broadening and blue shift. The blue shift is expected 
from the thermal expansion of the carbon-carbon bonds, 
which changes the tight-binding parameters [28|. On the 
other hand, the thermal broadening cannot be explained 
within the tight-binding model. To quantify the ther- 
mal broadening, we first fit the spectra within a single- 
particle model using the joint density of states for inter- 
band transitions, obtained from the SWM model, with 
T-dependent Lorentzian broadening. Figure [Hd) shows 
the extracted Lorentzian full width at half maximum 
(FWHM) r as a function of T for 20 and 30 T. Apart 
from a small finite line-width at T = 0, Fo ~ 5 cm^^, 
possibly due to disorder, F linearly depends on T. 

Within the single-particle picture, T only appears in 
the Fermi-Dirac distribution function, but this is a neg- 
ligible effect since both the initial and final states of the 
Raman process are far away from iJp, which resides in 
the n = bands. For example, for the (1,1) transition 
at 30 T, the electron and hole bands are ~65 meV (or 
~^75G K) away from E^- Thus, we need to take into 
account the interactions of the photo-excited electron- 
hole (e-h) pairs with some low-energy modes that would 
significantly change when T changes from 4 to 300 K. 
Specifically, since the linear-T broadening in Fig.[2Jd) im- 



plies a Bose-Einstein distribution at an energy scale much 
smaller than k^T, we only consider bosonic excitations 
whose characteristic energies are ^ 100 K. Hence, we 
consider two types of low-energy modes: i) particle-hole 
(p-h) excitations across Ep in the n = bands [Fig.[31^a)] 
and ii) acoustic phonons (the lowest optical phonons be- 
ing the shear modes, seen at ~5-6 meV in graphite [30|). 
We find that interactions with i) explain the observed T- 
linear broadening while interaction with ii) is too weak. 
The magneto-electronic Raman scattering matrix was 
previously calculated for SLG 31| and BLG [33 and can 
be readily generalized to graphite in the presence of B: 



R = Aj^Hiky^ k,)-^_^{ky, k,) 



(1) 



where A is the scattering amplitude, ky (kz) are electron 
momenta in the a6-plane (along the c-axis) , ^J^ creates an 
electron in the n = 1,2, 3, 4, ... bands, and ^-„ creates a 
hole in the n =— 1,— 2, —3, —4, ... bands. Both electrons 
and holes are massive at the bottom of the bands at the 
K point, i.e., to„ ^ for all n's, similar to BLG, but 
there is electron-hole asymmetry, i.e., mi ^ rn_i. 

Figure[3l(a) depicts the basic ingredients involved in the 
electron-electron interaction process we consider here, to- 
gether with dispersions calculated via the SWM model 
for the n = 0^ and ±1 bands at 20 T. The two lowest- 
energy bands (n = ) cross E'p, and the carriers near Ep 
have approximately linear dispersions. In the (1,1) pro- 
cess, an e-h pair is created in the n = ±1 bands, which 
interact with, and are thereby dressed with, multiple p-h 
excitations in the n = 0^ bands near Ep. As T is raised, 
the thermal smearing of the Fermi edge leads to stronger 
interaction between the massive e-h pair and the massless 
p-h pairs, and the peak broadens. This type of shake-up 
process was theoretically studied in carbon nanotubes at 
K [H,!!!!: a Id van Hove singularity, (w - Al-i/2 
predicted to become (w — A)~^/^+^" with a 



A)--'-, IS 
-^0.1 once 



the shake-up process is taken into account. 

We describe the n = 0~ electrons as_a Tomonaga- 
Luttinger liquid with the Hamiltonian [ij- 



m 



v-p / dz 



iptidzil^-R - il^Udzilj-L 



(2) 



where wp is the Fermi velocity and V'RfLi creates a particle 
near the right (left) Fermi point. The n = 0+ band can be 
described by a similar Hamiltonian, but with a different 
v-p. By approximating the energy dispersion near Ep as 
E oc k^, we rewrite Eq. ^ via bosonization: 



Ho = ^ldz[iy^r + {y9r] 



(3) 



where V0 = — 27r[(OR, + pi], V9 = 27r[pR — pl], and 
PR (pl) is the density operator for right-moving (left- 
moving) electrons. 
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FIG. 3: (Color online) (a) Electron-electron interaction pro- 
cess influencing the electronic Raman line-shape of graphite 
in high B. The (1,1) electron-hole pairs interact with, or shake 
up, the id Fermi sea in the lowest-energy Landau subbands, 
creating particle-hole pairs across the Fermi energy, (b) T 
dependence of the line-shape for the (1,1) transition at 20 T, 
together with fits (dashed lines) based on the model in (a). 



We assume that the photogenerated electrons (n = 1) 
and holes (n = — 1) interact with the n = Q~ conduction 
electrons separately. For the n = \ band, where electrons 
are massive, we can treat the electrons through: 



Hi 



rfz*i 



2m 



Ai 



*i, 



(4) 



where Ai is the band edge frequency and "^\ (^i) is 
the creation (annihilation) operator for the n = 1 band. 
We also assume that the interaction Hamiltonian only 
involves the total charge density, thus neglecting any 
backscattering and Umklapp scattering: 



Hint ^ — I dz 



ZTT 



(5) 



We write the effective Hamiltonian for the system as the 
sum of Eqs. dS])-©: H = H^ + Hi + Hint- Its diagonal- 
ization is a unitary transformation U'^ HU and has been 
previously solved by many authors |33l436| : 



C/t 



exp 



il- j dye{y)^\{y)-^i{y) 



(6) 



Under this transformation, the original interacting sys- 
tem can be mapped to a non-interacting one, and the 
massive-electron operator acquires an additional string 
operator, ^1(2:) — exp[— i7"'"6'(z)/7r]^i(z), where ^\ cre- 
ates a free electron in the n = 1 band. The massive n 
= 1 electron then gets dressed by the additional string 
operator, i.e., the n = 0~ conduction electrons adiabat- 
ically adjust to the massive electrons. Similarly, we can 
obtain a dressed expression for the massive hole. 

The spectral function can be obtained by calculating 
the imaginary part of the retarded Green's function (5|: 

G^(z,t) = -i0(t)([vl/t_^(z,i)*i(z,t),vl/|(O,O)vI/_i(O,O)]). 

(7) 
At zero T, Eq. ([7]) can be evaluated directly in real space. 
However, at finite T, one has to follow a different route. 
As the Green's function for the massive electron/hole and 
that for the conduction electrons are both straightfor- 
ward to obtain, the total Green's function can be written 
as a convolution of three Green's functions: 






-i9{t)[-iG<i{-z,-t)] 
x[iGf{z,t)]F{z,t), 

x5{0 - pi - p2)S{uj - wi - UJ2), 
f dpi f duji 

xGfi(pi -p,uji - w). 



where 



Fiz,t) 



(8) 



(9) 



(cxp[-i76'(x, t)] cxp[i76'(0, 0) 
We can express the spectral function in a universal form: 

'uj/T 



A{uj) = AT 



2Q-0.5 



F 



An 



■,a 



(10) 



where 






F{z,t) -- 


00 C« 

= ^ y] B[n -1- a, 1 - a]B[m 

n=Om=0 




xRe 


(2i)2" 




Jz- ^{m + n + a) 



a,l — a] 



(11) 



In Eq. (|lip there are two dimensionless parameters: 
Lu/T and a. The first implies that the spectral width 



linearly depends on T for a fixed a. The second can be 
understood by studying the T = asymptotic behavior 
of Eq. (IllD, an d comparing it with the previous zero-T 
results 



34| . It then becomes clear that: 

e(cj - A) 



Aiu) 



ex 



(W - A)l/2-2a 

where Q is the Heaviside function, 
otubes, a was estimated to be^ 0.1 



(12) 



For metallic nan- 
Mlsl. To fit our 



experimental data with our model, we use the true band 
structure instead of a parabolic approximation, by fit- 
ting the tail up to ^ 0.2(7r/c) from the K point. Figl3Jb) 
shows how well our model fits the data, giving q;(20T) 
w 0.016, a(30T) « 0.026, and a(40T) « 0.048. 

We now consider the longitudinal acoustic (LA) 
phonons in graphite, which can also couple to the massive 
electrons and holes. The properties of acoustic phonons 
can be described by five elastic constants [37|: Cn = 
1109 GPa, Cee = 485 GPa, C33 = 38.7 GPa, C13 = 
GPa, and C4 4 = 5 GPa. Unlike the case of optical 
phonons 19|, l26|, |38|, l39| , coupling with acoustic pho nons 
vanishes at the F point since H^p ~ ^/q |25l . |40| . |41|. 
where q is the phonon wavenumber. We then evaluate 
the thermal broadening of Raman peaks by calculating 
the imaginary part of the self-energy: 



H, 



k' ,k,q 



g^Hq)^\{ky + qy, k, + q,)'i'iiky, k,){bl^+bg') 



rjKq sm f 



V2A?, 



2NMujg- 



h{q) 



2l\q^ sin^ I 



2 r7i 

l%q^sin^i 



(13) 



-(iW 



•e)/4 



where Ib = {h/eBY^^ is the magnetic length, 77 ^ 2, and 
K ~ 1/3 [26|. To first order, we estimate the scattering 
rate through Fermi's golden rule: 



W, 



2tt 



Y,\{f\H.M'HE^~Ef). (14) 

/ 

When the momentum transfer during the scattering pro- 
cess is small (i.e., vq ^ kT), the phase space for phonon 
modes are g^( ^^/^t + ^ ^ ^) ^ qT ^ 0, and when the 
momentum transfer is large, the overlap between initial 
and final states has a factor exp(— g^Z^). For B = 30 T, 
Ib ^ ^ nm, which is at least one order larger than the 
carbon-carbon bond length. Thus, the contribution to 
scattering drops exponentially as the phonon modes move 
away from the F point (or, equivalently. with increasing 
energy). The calculated momentum-dependent scatter- 
ing rate is then given by 



W{k,) = A' 



de 



Sp s\v? I 



h\q,e) 



3^0 



1 1 

"-'^+2^2 



(15) 



where 

q = 
M' = 



2mlBV^ ll^ cos T ^8m'0 + Vi/V,^cos^e 



n cos"= 

ifK^ m Knit A| V2t; 
167r M Z| 72 K ^ "" 



!4.1 X 10"^ cm"^ 



This value leads to W{k,^) - lO'^ cm"! at 30 T and 
200 K, too small to explain the observed broadening, 
which requires the scattering rate to be ~10 cm~^. 
There are two reasons for the small KA' value: one is 
m/M ~ 10~^, and the other Knit/^l- The latter, i.e., 
the magnetic length suppression, is a unique aspect of 
this work, made possible by a high B. 

In summary, we studied electronic Raman scattering 
in graphite in a strong magnetic field up to 45 T, ap- 
plied along the c-axis. We observe a series of spectral 
lines, ascribed to inter-Landau-subband transitions, and 
each line exhibits strongly T-dependent line-shape. We 
developed a microscopic model based on the Tomonaga- 
Luttinger theory, with which we show that electron- 
electron interactions can explain the observed results, 
through the 'shake-up' process known in the problem of 
X-ray (or Fermi-edge) singularities. Specifically, we show 
that electron-electron interactions modify the van Hove 
singularity to the form {uj — A)^"~^/^ at T = 0. Our 
model accurately reproduces the observed T-dependent 
line-shape, determining a to be 0.016, 0.026, and 0.048, 
at 20, 30, and 40 T, respectively. 
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